arXiv:1502.01585vl [astro-ph.CO] 5 Feb 2015 


Astronomy & Astrophysics manuscript no. AOT LFFmapmaking © 2015 

February 6, 2015 


Planck 2015 results. VI. LFI mapmaking 

Planck Collaboration: R A. R. Ade^^, N. Aghanim^^, M. Ashdown^^’^, J. AumonF^, C. Baccigalupi^^ A. J. Banday^^’^, R. B. Barreiro^^, 

N. Bartolo^^’^^, E. Battaner^^’^^, K. Benabed^^’^^, A. BenoiF^, A. Benoit-Levy^^’^^’^^, J.-R Bernard^^’^, M. Bersanelli^^'^^, R Bielewicz^^’^’^^ 

A. Bonaldi^^, L. Bonavera^^, J. R. Bond^, J. BorrilF^’^^, R R. BoucheF^’^"^, M. Bucher^ C. Burigana"^^’^^"^^, R. C. Butler"^^, E. Calabrese^^, 
J.-E. Cardoso^^’^’^^, A. Catalano^^’^^, A. Chamballu^^’^^’^^, R.-R. Chary^"^, R R. Christensen^^’^"^, S. Colombi^^’^^, L. R L. Colombo^^’^"^, 

B. R CrilF"^’^, A. Curto^’^^, E. Cuttaia"^^, L. Danese^^ R. D. Davies^^, R. J. Davis^^, R de Bernardis^^, A. de Rosa"^^, G. de Zotti'^^’^^ 

J. Delabrouille^ C. Dickinson^^, J. M. Diego^^, H. Dole^^’^^, S. Donzelli"^^, O. Dore^"^’^, M. Douspis^^, A. DucouF^’^^, X. Dupac^^, 

G. Efstathiou^^, E. Elsner^^’^^’^^, T. A. EnBlin^^, H. K. Eriksen^^, J. Eergusson^^, E. Einelli"^^"^^, O. EornF^’^, M. Erailis"^^, E. Eranceschi"^^, 

A. ErejseF^, S. Galeotta"^^, S. GallF^, K. Ganga^ M. Giard^^’^, Y. Giraud-Heraud^ E. Gjerlpw^^, J. Gonzalez-Nuevo^^’^^ K. M. Gorski^"^’^^, 

S. Gratton^^’^^, A. Gregorio^^ A. Gruppuso"^^, E. K. Hansen^^, D. Hanson^^’^"^’^, D. L. Harrison^^’^^, S. Henrot-Versille^^, D. Herranz^^, 

S. R. Hildebrandt^"^’^, E. Hivon^^’^^, M. Hobson^, W. A. Holmes^"^, A. Hornstrup^"^, W. HovesF^, K. M. Huffenberger^^, G. Hurier^^, A. H. Jaffe^^, 

T. R. Jaffe^^’^, M. Juvela^^, E. Keihanen*^^, R. Keskitalo^^ K. KiiverF^ '^^, T. S. Kisner^"^, J. Knoche^^, M. Kunz^^’^^’^, H. Kurki-Suonio^^"^^, 

A. Lahteenmaki^'^^, J.-M. Lamarre^^, A. Lasenby^’^^, M. LattanzF^, C. R. Lawrence^"^, J. R Leahy^^, R. LeonardF^, J. Lesgourgues^^’^^’^^, 

E. Levrier^^, M. LiguorF^’^^, R B. Lilje^^, M. Linden-Vprnle^"^, V. Lindholm^^"^^, M. Lopez-Caniego^^’^^, R M. Lubin^^, J. E. Macias-Rerez^^, 

G. Maggio"^^, D. Maino^^"^^, N. MandolesF^’^^, A. MangillF^’^^, R G. Martin^, E. Martmez-Gonzalez^^, S. MasF^, S. Matarrese^^’^^’^^, 

R Mazzotta^^, R McGehee^"^, R R. Meinhold^^, A. MelchiorrF^ "^^, L. Mendes^^, A. Mennella^^ "^^, M. Migliaccio^^’^^, S. Mitra^^’^"^, L. Montier^^’^, 
G. Morgante"^^, D. Mortlock^^, A. Moss^^, D. Munshi^^, J. A. Murphy^^, R Naselsky^^’^"^, E. NatF"^, R NatolF^"^"^^, C. B. Netterfield^^, 

H. U. Nprgaard-Nielsen^"^, D. Novikov^^, L Novikov^^’^^, R Raci^F L. Ragano^^"^^, D. RaolettF^"^^, B. Rartridge'^F R Rasian"^^, G. RatanchonF 
T. J. Rearson^’^"^, O. Rerdereau^^, L. Rerotto^^, R Rerrotta^F V. Rettorino"^^, E. RierpaolF^, D. Rietrobon^"^, E. Rointecouteau^^’^, G. Rolenta"^"^"^, 

G. W. PratF^, G. Prezeau^’^"^, S. PruneF^’^^, J.-L. PugeF^, J. R Rachen^^’^^, R. Rebolo^^’^^’^^, M. Reinecke^^, M. Remazeilles^^’^^ F A. Renzi^^’^^, 
G. Rocha^"^’^, C. RossetF M. Rossetti^^"^^, G. Roudier^’^^’^"^, J. A. Rubino-Martm^^’^^, B. Rusholme^"^, M. SandrF^, D. Santos^^, 

M. Savelainen^^’"^^, D. Scott^^, M. D. SeifferF"^’^, E. R S. Shellard^^, L. D. Spencer^^, V. Stolyarov^’^^’^^, R. StomporF D. Sutton^^’^^, 

A.-S. Suur-Uski^^’"^^, J.-R SygneF^, J. A. Tauber^^, L. TerenzF^ "^^, L. Toffolatti^^’^^"^^, M. TomasF^"^^, M. Tristram^^, M. Tucci^^, J. Tuovinen^, 
L. Valenziano^F J. Valiviita^^’^F B. Van TenFF T. Vassallo^F R Vielva^F R. Villa^F L. A. Wade^F B. D. WandelF^’^^^F R. Watson^F 

I. K. Wehus^F D. Yvon^F A. ZaccheFF and A. Zonca^^ 

(Affiliations can be found after the references) 

Rreprint online version: Rebruary 6, 2015 


ABSTRACT 

This paper describes the mapmaking procedure applied to Planck LRI (Low Frequency Instrument) data. The mapmaking step takes as input the 
calibrated timelines and pointing information. The main products are sky maps of /, g, and U Stokes components. For the first time, we present 
polarization maps at LFI frequencies. The mapmaking algorithm is based on a destriping technique, enhanced with a noise prior. The Galactic 
region is masked to reduce errors arising from bandpass mismatch and high signal gradients. We apply horn-uniform radiometer weights to reduce 
effects of beam shape mismatch. The algorithm is the same as used for the 2013 release, apart from small changes in parameter settings. We 
validate the procedure through simulations. Special emphasis is put on the control of systematics, which is particularly important for accurate 
polarization analysis. We also produce low-resolution versions of the maps, and corresponding noise covariance matrices. These serve as input in 
later analysis steps and parameter estimation. The noise covariance matrices are validated through noise Monte Carlo simulations. The residual 
noise in the map products is characterized through analysis of half-ring maps, noise covariance matrices, and simulations. 

Key words, methods: data analysis - cosmic microwave background - numerical methods 


1. Introduction 

This paper is one of a set associated with the 2015 release of 
data from the Planck^ mission (Planck Collaboration I 2014). It 
describes the mapmaking procedure applied to time-ordered data 


* Corresponding author: E. Keihanen elina.keihanen® 
helsinki.fi 

^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from Erance and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 


from the Planck Low Frequency Instrument (LFI). Mapmaking 
is the next step in the LFI pipeline after calibration and dipole 
removal. It is followed by bandpass correction and component 
separation. A description of the complete data analysis pipeline 
is given in Planck Collaboration II (2015). 

The mapmaking procedure applied to 2013 data was dis¬ 
cussed in Planck Collaboration II (2014). In this work we present 
a thorough validation on the procedure, and set focus on the 
control of systematics (Planck Collaboration III 2015). For the 
first time we present polarisation results. The mapmaking ap¬ 
proach adopted for HFI is described in Planck Collaboration 
VIII (2015). 

Mapmaking takes as its input the calibrated timelines, from 
which the cosmological and orbital dipole signals have been re- 
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moved. Also, an estimate for the Galactic straylight is subtracted 
from the timelines prior to mapmaking, since this is difficult to 
correct for at map level. Calibration and removal of stray light 
are described in Planck Collaboration V (2015). 

In addition to the time-ordered data (TOI), required inputs 
are, the radiometer pointing (6, 0,0) for each sample, informa¬ 
tion of the polarization sensitivity of each radiometer, and noise 
description. Output consists of sky maps of temperature and 
Stokes Q and U polarization, and a description of residual noise 
in them. Effects of beam shape are not corrected for at the map¬ 
making level. The LFI beams and the associated window func¬ 
tions are described in Planck Collaboration IV (2015). 

The time-ordered data are contaminated by correlated 1 // 
noise. The aim of the mapmaking step is to remove the correlated 
noise as accurately as possible, while simultaneously minimis¬ 
ing systematic errors. Various mapmaking algorithms in the con¬ 
text of Planck are compared in Poutanen et al. (2006); Ashdown 
et al. (2007a,b, 2009). 

The mapmaking step is computationally challenging, due to 
the huge amount of data. For optimal removal of noise, all of the 
time-ordered data must be held in computer memory at the same 
time. This leads to a large memory requirement (of the order of 
a TB for 70 GHz). The size of the output map is much less (of 
the order of 100 MB per map). Mapmaking is thus also a data 
compression step. 

Providing a statistical description of the residual noise 
present in the maps, in the form of pixel-pixel noise covari¬ 
ance matrices, is part of the mapmaking step. Covariance matri¬ 
ces are used in a number of different places in the data-analysis 
pipeline. These matrices are impossible to employ at the native 
map resolution because of resource limitations. At present, the 
low-resolution products can only be used efficiently at resolution 
A^side = 16 or lower. Consequently, we also produce correspond¬ 
ing low-resolution maps from the high-resolution maps. 

This paper is organized as follows. In Sect. 2 we provide 
an overview of the mapmaking procedure, with emphasis on the 
control of systematics. Validation of the procedure with simula¬ 
tions is presented in Sect. 3. Production of low-resolution maps 
and related noise covariance matrices are discussed in Sect. 4. 
In Sect. 5 we give an overview of the mapmaking products and 
characterize the residual noise in them. A summary and conclu¬ 
sions are given in Sect. 6. 

2. Mapmaking procedure 

2.1. Madam 

LFI maps were produced by the Madam mapmaking code, ver¬ 
sion number 3.7.4. The algorithm is based on the destriping 
technique, enhanced with a noise prior. The destriping tech¬ 
nique as a tool for mapmaking is discussed in Maino et al. 
(2002), Keihanen et al. (2004), and Kurki-Suonio et al. (2009). A 
noise prior allows us to extend the destriping technique to very 
short baseline lengths, allowing more accurate noise removal. 
The Madam algorithm is described in Keihanen et al. (2005) and 
Keihanen et al. (2010). The code is the same as used for the 2013 
release, but the parameter settings have been altered. The 2013 
mapmaking procedure was discussed in Planck Collaboration 
II (2014). Here we give an overview of the method for conve¬ 
nience. We focus on changes made after the 2013 release, and 
on aspects relating to polarization. 

Consider a time-ordered data stream 

y = s n' (1) 


where s represents the sky signal and n' is noise. In the destrip¬ 
ing technique, the correlated noise component is modelled by 
a sequence of offsets, here called baselines. The same offset ap¬ 
plies for N adjacent samples in time-ordered-data. In the conven¬ 
tional language N is referred to as the “baseline length,” while 
“baseline” alone refers to the value of the offset. 

The baselines are assumed to be the result of a random pro¬ 
cess. On top of the correlated noise, pure white noise with known 
variance is assumed. Formally we can write 

n' = Va + n. (2) 

Vector a represents the baselines, and F is formally a matrix that 
spreads the baselines into the time-ordered data. Vector n repre¬ 
sents the white noise component. The noise model ignores corre¬ 
lated noise at frequencies that exceed the inverse of the baseline 
length. The technique is suitable for 1 //-type noise that dies out 
at high frequencies. The baseline length should be chosen to be 
short enough to capture all frequencies below or comparable to 
the knee frequency of the 1// spectrum. 

The signal component is modelled as 

s = Pm, (3) 

where m is a pixelized sky map, and P is a pointing matrix. Map 
m consists of three components, representing temperature and Q 
and U polarization. We can thus write for sample /: 

Si = T{pi) + Q(pi) cos(2i/r,) + U(pi) sin(2i/r,). (4) 

Here pi is the pixel hit by sample i, and if/i is the rotation an¬ 
gle between the local meridian and the vector that defines the 
direction of polarization sensitivity of the beam. 

Each sample is assigned entirely to the pixel in which the 
beam centre falls. Beam shape is not taken into account in this 
context. We are also ignoring cross-polarization, which for LFI 
is below 1% (Planck Collaboration III 2015). Formally then, P 
is a sparse matrix with three non-zero elements in each row. The 
factors cos(20/) and sin(20/) are included in P. Limitations of 
the model and their impact on mapmaking are discussed in Sect. 
2.3. 

In the conventional destriping technique no prior information 
on the baselines is used. On the other hand, if the noise spectrum 
is known, one can construct a noise prior 

Ca = {aa^). (5) 

The noise prior provides the extra constraint that allows us to 
extend the destriping technique to very short baselines. We con¬ 
struct the noise prior from the known noise parameters (knee fre¬ 
quency, white noise variance, and spectral slope), as described 
in Keihanen et al. (2010). It can be shown that the destriping 
solution for the map gives the maximum-likelihood map in the 
limiting case where the baseline length approaches one sample 
(Keihanen et al. 2005). 

With the assumptions listed above, the baseline vector a can 
be solved from the linear system of equations 

(FX‘ZF + C:i)a = FTc;'Z3^, (6) 

where 

z = \-p(p"^c-^pr^p^c-\ ( 7 ) 

Here Cw is the white noise covariance, formally a diagonal ma¬ 
trix with the white noise variance on the diagonal. The final de- 
striped map is then constructed as 

m = (P'fc;‘pr‘p'fc;‘0'-Fa), (8) 
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where a is the baseline solution from (6). We use the hat symbol 
to indicate that a and m are estimates for the true a and m. The 
“binned map” is constructed from the timeline the same way 
without baseline removal 

mb = (P^C;‘P)-‘pTC-y (9) 

We include the same set of radiometers and the same time 
range in destriping phase and in map-binning phase. In this re¬ 
spect our approach is different from the one adopted for HFI 
(Planck Collaboration VIII 2015). Our approach has the benefit 
that single-year or single-survey maps, or maps constructed from 
different radiometer sets, have independent noise contributions. 
The price to pay is that the baseline solution in the partial maps 
is less accurate than in the full survey maps. 

Adding a uniform offset to the time-ordered data does not 
have an impact on the baseline solution. This can be seen by 
inserting a unity vector j = 1 to the right-hand side of Eq. 
(6). A unity vector can be expressed as 1 = Pnti, where m\ is 
an (/, 2, U) map with temperature component equal to one and 
with vanishing polarization. The operation ZPnti yields zero, 
and consequently to a zero baseline solution. Based on similar 
arguments, destriping cannot determine the absolute monopole 
of the temperature map, since this is indistinguishable from an 
offset in the timeline. 

We do not impose an explicit constraint on the sequence of 
baselines to avoid the degeneracy between the monopole of the 
maps and an offset in TOI. Instead, we make use of the fact the 
conjugate gradient iteration method works for singular systems 
as well, under the condition that the same eigenmodes vanish on 
both sides of the equation. The raw output temperature map will 
in general have an arbitrary non-zero mean that depends on the 
input signal and noise in a non-trivial way. A monopole is sub¬ 
tracted from each released temperature map in post-processing, 
as described in Planck Collaboration II (2015). 

Similar arguments do not apply to an offset difference be¬ 
tween radiometers. In contrast to a common offset, a differential 
offset cannot be presented as a map, because it does not obey 
the rotation properties of the 2, U maps. A differential offset is 
distinguished from a map by the destriping procedure, and does 
not leak into polarization, as long as there are enough crossing 
points where the same detector scans the same pixel in different 
orientations. 

The destriping technique produces the final map through a 
procedure in which one first solves for the baselines, and then 
bins the map from the data stream from which the baselines have 
been removed. This two-step procedure has a benefit over direct 
maximum-likelihood methods, because it enables us to use tech¬ 
niques to make the map diverge from the maximum-likelihood 
solution and thus have a better control of systematics. For in¬ 
stance, one can apply a mask in the destriping phase while still 
binning the final map to cover the whole sky, in order to reduce 
signal error due to strong signal gradients. We discuss this in 
more detail in Sect. 2.3. 

2.1.1. Baseline length and resolution 

The baseline length is a key parameter in the Madam mapmak¬ 
ing algorithm. The chosen baseline length is a trade-off between 
computational burden and optimal noise removal. We have se¬ 
lected baselines that are well below the timescale given by the 
knee frequencies of each channel. In the case of 44 GHz and 
70 GHz we have chosen to use one second long baselines. For 
30 GHz, where the knee frequencies are higher, we have used 


shorter baselines of 0.25 s. This is different from the 2013 re¬ 
lease, where we used one second baselines for all channels. 
The exact baseline lengths and channel sampling frequencies are 
given in Table 1 . The FFI knee frequencies are given in Table 2 
along with other noise parameters. 

Madam maps follow the HEALPix^ pixelization scheme, 
where map resolution is defined by parameter Aside- The total 
number of pixels in map is 12^^^.^^ (Gorski et al. 2005). Map 
resolution was chosen according to the beam width. For the ma¬ 
jority of maps we used resolution Aside=1024. At this resolution, 
the average pixel size corresponds to 1/10 of the beam FWHM 
at 30 GHz or 1/4 at 70 GHz. A subset of 70 GHz maps was re¬ 
done with resolution Aside=2048. All computations were carried 
out in Galactic coordinate frame. Pixel sizes are compared with 
beam widths in Table 2. 

2.1.2. Missing data 

Periods of stable pointing are regularly interrupted by manoeu¬ 
vring periods, where the telescope is repointed according to the 
scanning strategy. These repointing periods are discarded from 
mapmaking. In addition there are other data sections that must 
be discarded, for instance due to missing data, gain saturation or 
crossing of a bright object (Planck Collaboration II 2015). We 
accomplish this by formally setting the white noise variance Cw 
to infinity for the flagged samples. The procedure is different 
from removing the flagged data sections and appending the re¬ 
maining sections end-to-end. As a result, the flagged samples do 
not contribute to the output map, but the correct noise correlation 
is preserved across the gaps. 

2.1.3. Single-survey maps and handling of missing pixels 

The 2015 release includes polarization maps for individual sur¬ 
veys (6 months of data), which do not cover the whole sky. 
Obviously, pixels that are not observed at all, are not included 
in the maps. In addition there are pixels that are scanned by a 
single horn. The polarization analysis of those pixels requires 
special attention. To reliably solve the /, Q, and U components 
for a pixel, the pixel has to be scanned by two horns with or¬ 
thogonal polarization sensitivity, or, alternatively, twice by same 
horn with sufficiently different orientations of the beam. If this 
is not the case, we discard the pixel from the map. Even these 
pixels, however, can be used for destriping as we now describe. 

We assess the goodness of a pixel through the matrix 

P^CP. (10) 

This is a block-diagonal matrix with 3x3 blocks along the diag¬ 
onal. Each 3x3 block corresponds to one pixel, and represents 
the inverse of the white noise covariance within the pixel. The 
matrix gives an accurate diagnosis of a pixel’s usefulness for po¬ 
larization analysis. Pixels with insufficient polarization coverage 
yield a singular or badly scaled matrix. The inverse of the pixel 
matrix appears in the destriping Eqs. (6) and (7). It is obvious 
that a singular matrix leads to a non-stable solution for the map. 
However, in Eq. (6) the pixel matrix appears sandwiched be¬ 
tween two occurrences of the same pointing matrix, which lack 
the same eigenmodes. Matrix Z is thus finite even if the pixel 
matrix is singular. This can also be seen from the fact that Z is a 
projection matrix, which means that all its eigenvalues are 1 or 
0 . 

^ http://healpix.sourceforge.net 
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Table 1. Frequency-specific mapmaking parameters and related information. From left: sampling frequency; knee frequency range; 
and chosen baseline length in seconds and as a number of samples. The three rightmost columns show the FWHM range of the main 
beams of the radiometers of the channel, the chosen destriping resolution as the A^side parameter of HEALPix, and the average size 
of one pixel. 


Baseline length Resolution 


Channel 

/amp [Hz] 

An [Hz] 

s 

samples 

FWHM [arcmin] 

A^side 

arcmin 

30 GHz 

32.508 

0.043-0.175 

0.246 

8 

33.0-33.2 

1024 

3.44 

44 GHz 

46.545 

0.020-0.088 

0.988 

46 

23.0-30.8 

1024 

3.44 

70 GHz 

78.769 

0.006-0.059 

1.000 

79 

12.8-13.5 

1024/2048 

3.44/1.72 


If the pixel matrix is singular, but not identically zero, we 
invert the non-singular eigenmodes. That way we can also uti¬ 
lize pixels that have been observed by a single horn. In the 
final map-binning phase we cannot use this same procedure. 
Instead, the badly determined pixels are excluded. We use as re¬ 
jection criterion the reciprocal condition “rcond” value of the 
pixel matrix. For the 2015 release we chose the conservative 
limit rcond > 0.01. Statistics on the sky coverage of LFI maps is 
provided in Sect. 5. 

2.1.4. Single-horn and single-radiometer maps 

In order to decompose the signal into /, Q, U components, it is 
necessary to have at least three measurements from the same 
pixel, in different polarization angles. In a limited region around 
the ecliptic poles it is possible to determine the polarization sig¬ 
nal already from one LFI horn, since the region is scanned in 
different angles by the same radiometer pair, but for most of the 
sky we need observations from two horns (four radiometers). 

We have produced temperature maps from each horn’s or 
radiometer’s data. The same formulation as above applies, but 
the pointing matrix P only contains the temperature column. 
Consequently, the matrix in Eq. (10) becomes diagonal and is 
trivial to invert. 

In the case of horn maps, the polarization signal largely can¬ 
cels out when the two radiometers are combined, and the map 
gives a reliable estimate of the temperature signal. The single¬ 
radiometer maps, instead, are a mixture of temperature and po¬ 
larization signals, which must be kept in mind when analysing 
these maps. 


2.2. Noise prior 

We use a prior C“^ for the noise baselines to constrain the de¬ 
striping solution. Matrix Ca gives the expected correlation be¬ 
tween noise baselines, as in Eq. (5). The details of the compu¬ 
tation from a given input spectrum are given in Keihanen et al. 
(2010). The prior makes it possible to use very short baselines, 
and thus to model the noise accurately. 

We construct the prior from the noise parameters (cr, 
which are listed in Table 2. The noise spectrum is given by 


P(f) = 


— 

r\r ’ 


( 11 ) 


where /samp is the sampling frequency (see Table 1). Below a 
cut-off frequency /min =1/(3600 s), the spectrum is assigned a 
constant value P(fmin), to avoid numerical problems related to 
the divergence at / = 0. 


On top of the correlated noise component there is assumed 
to be pure white noise with constant spectrum 

PUf) = ( 12 ) 

/samp 

The noise prior is based on the correlated component of Eq. (11) 
only, since this is the component modelled by the baselines. 

We assume that the noise properties do not change signifi¬ 
cantly with time, and apply the same prior throughout the mis¬ 
sion. We also assume that noise is uncorrelated between ra¬ 
diometers, and from one pointing period to another. Under these 
assumptions the correlation matrix Ca has a simple form and can 
be applied efficiently using Eourier techniques. 

2.3. Handling of systematic effects 

Signal variations that are not captured by the signal model (Eq. 
(3)) give rise to additional uncertainty in baseline determination, 
referred to as “signal error”. There are three main sources for 
signal error: signal variations within a pixel; differences in ra¬ 
diometer frequency response (bandpass mismatch) (Leahy et al. 
2010); and beam shape mismatch. Cross-polarization in this con¬ 
text can be understood as part of beam mismatch. 

Other sources include variable sources and gain estimation 
error from the preceding calibration step. All these have the 
effect that the same sky pixel, when observed by different ra¬ 
diometers or in different orientation of the beam, or just with the 
beam centre in slightly different position, yields different signals. 
Madam interprets this as additional noise and fits baselines to it. 

Note that a symmetric beam component, with equal width at 
all radiometers that participate in the mapmaking process, does 
not give rise to signal error, since it does not create inconsis¬ 
tencies between the measurements. Instead, it the beam smooth¬ 
ing is just interpreted as a property of the sky, and can be cor¬ 
rected for in a power spectrum analysis by applying an appropri¬ 
ate beam window function. 

The same error sources also cause leakage from the temper¬ 
ature signal into the weaker polarization maps. The relative im¬ 
portance of systematic effects, and their effect on LEI maps, is 
discussed in Planck Collaboration III (2015). 

The effects of beam shape mismatch and cross-polarization 
are not taken into account in the mapmaking step, but they are 
handled in the later data processing steps through application of 
a beam window function (Planck Collaboration IV 2015). Other 
methods have been also developed for full beam deconvolution 
(Keihanen & Reinecke 2012). 

Ideally, optimal noise removal is achieved by destriping the 
data according to the maximum-likelihood solution of Eqs. (6- 
7), with Cw equal to the actual white noise covariance, and using 
all available data. However, this does not take into account the 
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Table 2. Noise parameters for LFI radiometers. Values for /knee,y^, and cr given here were used for the construction of the noise prior. 
We give also radiometer weights for pure noise weighting and for the horn-uniform weighting scheme. Horn-uniform weighting 
was used for actual mapmaking. 



Knee Frequency Anee [mHz] 

Slope p 

White noise 

o-[mK] 

Weight 

^ (norm.) 

Horn-unieorm 


M 

S 

M 

S 

M 

S 

M 

S 

M&S 

70 GHz 

LFI-18 

14.82 

17.78 

-1.060 

-1.180 

4.553 

4.146 

0.1676 

0.2226 

0.1945 

LFI-19 

11.72 

13.71 

-1.207 

-1.111 

5.144 

4.926 

0.1313 

0.1577 

0.1454 

LFI-20 

7.96 

5.67 

-1.198 

-1.298 

5.212 

5.507 

0.1279 

0.1262 

0.1283 

LFI-21 

37.89 

13.27 

-1.247 

-1.205 

4.003 

4.971 

0.2168 

0.1548 

0.1811 

LFI-22 

9.68 

14.80 

-1.424 

-1.237 

4.356 

4.715 

0.1831 

0.1721 

0.1790 

LFI-23 

29.65 

59.03 

-1.073 

-1.211 

4.476 

4.790 

0.1734 

0.1667 

0.1716 

44 GHzGHz 

LFI-24 

26.78 

88.30 

-0.942 

-0.908 

3.159 

2.734 

0.3163 

0.3421 

0.3279 

LFI-25 

20.07 

46.37 

-0.845 

-0.904 

2.834 

2.698 

0.3931 

0.3513 

0.3739 

LFI-26 

64.42 

68.19 

-0.918 

-0.758 

3.295 

2.887 

0.2907 

0.3066 

0.2982 

30 GHz 

LFI-27 

174.53 

108.79 

-0.927 

-0.907 

1.605 

1.729 

0.5602 

0.4716 

0.5167 

LFI-28 

130.08 

43.08 

-0.931 

-0.900 

1.812 

1.633 

0.4398 

0.5284 

0.4833 


systematic effects, which are particularly important for accurate 
polarization analysis. 

The two-step nature of the destriping procedure enables us to 
diverge from the ML solution to reduce systematics, at the cost 
of slightly increased residual noise. Residual noise is easier to 
account for, and is thus less harmful in subsequent analysis, than 
systematic effects. 

2.3.1. Destriping resolution 

Signal error arising from high signal gradients, and consequent 
variations inside a pixel, can be reduced with an increased de¬ 
striping resolution. This leads to slightly higher residual noise, 
as there are more unknowns to be solved for in the same amount 
of data. 

The destriping resolution can be chosen independently from 
the final map resolution. In high-resolution mapmaking we have 
used the same resolution (Vside=1024 or Vside=2048) both for 
destriping and for map binning. The choice of resolution plays a 
more important role in the construction of low-resolution maps 
(Sect. 4). 

2.3.2. Destriping mask 

Signal error originates mostly in the Galactic plane where signal 
gradients as well as bandpass effects are largest. We can greatly 
reduce the signal error by applying a Galactic mask in the de¬ 
striping phase. This is done by formally setting the white noise 
covariance to infinity for the samples under the mask, thus en¬ 
suring consistent treatment throughout. All samples are again 
included when binning the final map. Using a mask leaves fewer 
crossing points between scanning rings to be used for baseline 
determination. The residual noise is therefore expected to in¬ 
crease slightly. 

The choice of the destriping mask is a trade-off between min¬ 
imization of signal error and keeping enough data for accurate 
destriping. We use the same masks that were used in Planck 
Collaboration V (2015). The masks are illustrated in figure 1. 


The effect of masking in the destriping phase on the signal error 
is illustrated in Fig. 2. 

2.3.3. Radiometer weighting 

Individual radiometers are weighted according to Cw when data 
from different radiometers ar combined into one map according 
to Eqs. (6-8). In the maximum-likelihood solution Cw equals the 
white noise variance of the time-ordered data. We call this the 
“noise weighting”. 

The purpose of “hom-uniform weighting” is to reduce leak¬ 
age from temperature to polarization. We replace the white noise 
covariance Cw for a radiometer by the average of the covariances 
of the two radiometers of a horn. The weights that go on the di¬ 
agonal of Cw are constructed as 


where ctm and cts are the white noise rms values for the two 
radiometers sharing a horn. The same weight is applied to both 
radiometers. We are neglecting variations in white noise level. 
C“^ can be regarded as a radiometer-specific weighting factor. 
The values of C“^ for both horn-uniform-weighting and pure 
noise weighting are given in Table 2. The numbers are based 
on the (T values listed in the same table. Since the same matrix 
Cw appears on both sides in the destriping equations, the units 
and overall normalization are irrelevant. The values given in the 
table have been normalized so that the weights of all radiometers 
of the frequency channel add up 1. 

At the same time we make the flags uniform, so that if a 
sample is flagged as unusable at one radiometer of the horn, we 
flag the corresponding sample for the other radiometer as well. 
This, together with the horn-uniform weighting, ensures that the 
polarization map is solved solely from the difference of the two 
timestreams. This has the advantage over noise weighting that 
any instrumental effects that are equal for both radiometers effec¬ 
tively cancel out. In particular this is expected to apply to beam 
effects, since the beam shapes are similar for two radiometers 
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Fig. 1. Destriping masks for 30 GHz (top), 44 GHz, and 70 GHz 
(bottom). 


sharing a horn (Planck Collaboration IV 2015). The cancella¬ 
tion is not perfect, however, because the polarization angles of 
the two radiometers are not exactly orthogonal, and because the 
noise priors are different. 

The effect of each of these procedures on signal error and 
on residual noise are assessed through simulations, which are 
described in Sect. 3 

2.4. White noise covariance 

When noise weighting is applied, Cw equals the actual white 
noise covariance, and the inverse of Eq. (10) directly gives the 
covariance of residual white noise in a pixel. The diagonal el¬ 
ements give the variance of /, Q, U components, and the off- 
diagonal elements represent cross-correlations between the com¬ 
ponents. 

The horn-uniform radiometer weighting complicates the pic¬ 
ture. If Cn is the actual white noise variance in the time-ordered- 
data, and Cw the one used for weighting, the residual white noise 
in the map domain is obtained as 

Cwn^ap = (P^c;'P)-‘pTc;'c„c;‘P(pTc;‘P)-‘. (i4) 

Obviously, this reduces to the inverse of the expression in Eq. 
(10) when Cw = Cn- It can be shown that Eq. (14) obtains 
its minimum for fixed Cn when Cw = Cn, and this minimum 



Fig. 2. Effect of destriping mask on 70 GHz signal error. We 
show the simulated signal error in temperature, without (upper) 
and with (lower) applying a Galactic mask. 


value is (P^C“^P)“^ Noise weighting thus minimizes the resid¬ 
ual white noise. We compute the white noise covariance (14) for 
every map produced, and provide it as auxiliary information. 


2.5. Half-ring maps 


Half-ring maps are used for chacterization of total residual noise. 
We split each pointing period in two, and use the first halves for 
the first half-ring map, and the second halves for the second half¬ 
ring map. Very long pointing periods are first split into sections 
of one hour at maximum. These are further split into two half¬ 
ring parts. We then construct maps from the two halves with the 
same procedure as for the full maps. 

The two half-ring maps have nearly identical signal compo¬ 
nent, but independent noise contributions, apart from the very 
lowest frequency components. We take the difference between 
the two maps to cancel the signal. A correction is applied to 
account for the small differences in the hit count distribution be¬ 
tween the two halves. The “half-ring noise map” in a given pixel 
is constructed as 


nih 


mh\ - mi 

Wh 


(15) 


where mhi and mh 2 are the two half-ring maps, and Wh is a 
weighting factor given by 


Wh 


(^hl + ^h2) 




1/2 


(16) 


Here rihi and are the hit counts in the pixel under considera¬ 
tion. If the hit distributions are identical, this reduces to Wh = 2. 

The half-ring noise map consists of noise with similar prop¬ 
erties as the residual noise in the full maps, apart from the very 
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longest time scales (over an hour) that are cancelled in the pro¬ 
cess. The noise map provides a valuable tool for the characteri¬ 
zation of residual noise in the maps, since it is constructed from 
the same data as the maps themselves, and is independent of any 
noise model. 

We expect a small error in the noise estimate based on the 
half-ring difference maps to come from the regions of strong 
temperature gradients. In these regions the signal can “leak” into 
the half-ring difference maps. We call this effect the gradient 
leakage. It is caused by the slightly uneven distribution of the 
hits within each pixel. Although in the single-survey half-ring 
maps these effects are visible by eye in the Galactic region, their 
effect on the noise spectrum derived from the map is negligible, 
because they only affect a few pixels with strongest signal gra¬ 
dients. In the full mission maps the gradient leakage is further 
mitigated by the much larger hit count. 

3. Validation of the mapmaking procedure 

3.1. Simulations 

We have performed a series of simulations to assess the effect 
of various parameter choices on map accuracy. Our goal here is 
not to build an error budget, but to justify the parameter choices 
done in the mapmaking process. The impact of various system¬ 
atic effects on the Planck data products is assessed in Planck 
Collaboration III (2015). 

The simulations include noise with realistic noise parame¬ 
ters, and CMB-rforeground with realistic beam shapes. In all 
cases we used realistic LFI radiometer pointings and flags, and 
constructed /, Q, U maps at HEALPix resolution Aside=1024. 

3.1.1. Noise simulations 

Madam includes an internal noise generator that we used for noise 
simulations. The internal generator has an advantage over exter¬ 
nal codes, because it avoids the step of writing the generated 
timeline on disk, thus reducing the computational cost of the 
simulation. 

The code generates white and 1 // noise according to a given 
noise spectrum, inserts it in place of the time-ordered data, and 
processes it into a map in the same way as actual data. The re¬ 
sulting map is the “residual noise map”. 

The residual noise is dominated by white noise. It therefore 
makes sense to look at the correlated residual noise component 
in isolation. The residual noise map can be split into two inde¬ 
pendent components 

= fflwn + fflcRN- (17) 

Here represents the binned white noise map, constructed as 

m,,„ = (pTC-‘P)-‘pTC-Vw« (18) 

where j^n is the separately generated white noise timeline. The 
binned white noise map does not depend on the destriping pro¬ 
cedure, apart from the radiometer weighting. 

We then And the “correlated residual noise” (CRN) map by 
subtracting from the residual noise map. The CRN compo¬ 
nent includes the correlated noise residual that remains after de¬ 
striping, and the error in baseline determination. It can be shown 
that as long as the destriping resolution is equal to the map res¬ 
olution or exceeds it, the two components are statistically inde¬ 
pendent, i.e., 

(/w^„mcRN> = 0. (19) 


This makes the rms of the CRN map a convenient figure of 
merit. It shows the effects of changes in the mapmaking pro¬ 
cedure more clearly than the full noise. The total residual noise 
rms is obtained as the quadrature sum of the CRN and white 
noise values. 

In Monte Carlo simulations we generate the white noise 
stream independently from the correlated 1 // noise stream. This 
makes it possible to construct the binned white noise map and to 
subtract it from the total noise map to obtain the CRN map. With 
flight data this distinction is not possible. 

When generating noise we used a time-dependent noise 
model, where the time span of the mission was split into 5 
sections, each of which had a different noise spectrum. At the 
time of constructing the simulations this was the most accurate 
model available. When destriping the data we used a single noise 
prior constructed from a constant set of noise parameters. This 
is equivalent to the procedure used for actual mapmaking. To re¬ 
duce scatter, we have generated 10 noise realizations for each 
combination of mapmaking parameters. The values shown are 
obtained as quadrature sum of the individual rms values. 

3.1.2. Signal simulations 

We perform signal simulations where we generate time-ordered 
data starting from a sky model. We used as input the FFP7 sky 
model (Planck Collaboration XII 2015) which was the version 
available at the time of our simulations. Differences with respect 
to the release version FFP8 are small, and unlikely to affect the 
results. Each radiometer has its own input sky, integrated with 
the radiometer’s frequency response. We use the tools of the 
Levels package (conviqt, multimod) (Reinecke et al. 2006), 
to convolve the sky with the beam, and to scan a time-ordered 
data stream from the input sky according to the real radiometer 
pointing. We include the main and intermediate beam compo¬ 
nents, but not the far sidelobe component, since we assume that 
the latter is removed in the calibration process. We apply flags, 
and run the data through the mapmaking procedure. 

The “signal error rms” is computed as the rms of the differ¬ 
ence between destriped and binned signal maps constructed from 
same data stream. It has become customary to call this quantity 
signal error, though its interpretation as pure error is question¬ 
able. For instance, destriping may remove some of the system¬ 
atic polarization signal that arises from bandpass mismatch, but 
this shows up as an increased signal error. The difficulty here is 
that the binned map itself is subject to bandpass and beam mis¬ 
match, and there is not a unique “true” map to which we could 
compare the destriped maps. Signal error is rather a measure of 
how much the destriping process modifies the signal, compared 
to naive binning. 

Signal error rms is not a useful figure of merit when we 
consider the effect of radiometer weighting, because the chosen 
weighting scheme affects the binned map as well. We performed 
another series of signal simulations where we included only un¬ 
polarized CMB in the input sky. All signal in the output polar¬ 
ization is thus leakage from temperature. Analysis of the polar¬ 
ization maps gives information on how the chosen radiometer 
weighting scheme impacts the leakage. 

3.2. Results 

3.2.1. Baseline length and mask 

In the first series of simulations we varied the baseline length 
and examined its effect on residual noise and signal error. The 
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Fig. 3. Correlated residual noise (CRN) as a function of baseline 
length for all LFI frequencies. CRN is computed as the rms of 
the difference between the destriped map and the binned white 
noise map. We plot CRN for 1 (orange), Q (green), and U (blue) 
as a function of baseline length. Solid and dashed lines show the 
result with and without the destriping mask, respectively. The 
circles indicate the baseline value used in actual mapmaking. 


results are plotted in Figs. 3 and 4. We show the CRN and signal 
error RMS for temperature and 2, U polarization. Because the 
amplitudes of residual noise and signal error are very different, 
we plot them in different panels. The total residual noise is a sum 
of the CRN component, the binned white noise component, and 
the signal error. Since the three components are independent, the 
square of total residual noise rms is obtained as a quadrature sum 
of the individual rms values. The rms of the binned white noise 
component, again averaged over 10 noise realizations, are found 
in Table 3. 

As expected, residual noise decreases with decreasing base¬ 
line width. The values used in actual mapmaking are indicated 
by circles. The 0.25 s baseline length used at 30 GHz appears to 
be a bit of overkill, since the residual noise level changes only 
very little below one second (0.2 yuK). 

At 30 GHz and 70 GHz, Q and U have nearly identical noise 
levels. This is a consequence of the focal plane design. LFI horns 


I r\ 



10 ° 10 ^ 10 ^ 


Baseline length [s] 

Fig. 4. Signal error as a function of baseline length for all LFI 
frequencies, for I (orange), Q (green), and U (blue). Signal error 
is computed as the rms of the difference between destriped and 
binned maps. Solid and dashed lines show the result with and 
without the destriping mask, respectively. 


are arranged in pairs that have their direction of polarization sen¬ 
sitivity at 45° angles to each other, so that when one horn is 
measuring Q, the other is measuring U, and both components 
are recovered with same accuracy. An exception from this is the 
44 GHz channel, where horn LFI-24 does not have a counter¬ 
part. This leads to a difference in noise levels between Q and U 
at 44 GHz. 

The level of signal error is a result of complicated inter¬ 
play between destriping, phenomena causing signal error, and 
instrument design. Different phenomena dominate at different 
frequencies. Signal variations within a pixel and beam shape 
mismatch contribute more strongly to the signal error in tem¬ 
perature than in polarization, because the beam patterns of two 
radiometers of a horn are similar, so that the effect partly cancels 
in polarization. The same does not apply to bandpass mismatch. 

We ran the same series of simulations with and without de¬ 
striping mask. As discussed in Sect. 2.3, the mask is expected to 
reduce the signal error arising from high signal gradients. The 
masks used are depicted in Fig. 1. Figure 2 illustrates the effect 
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Fig. 5. Correlated residual noise (CRN) as a function of destrip- 
ing resolution for all LFI frequencies. CRN is computed as the 
rms of the difference between the destriped map and the binned 
white noise map. We plot CRN for I (orange), Q (green), and U 
(blue) as a function of destriping resolution. The circles indicate 
the resolution (A^side=1024) used in actual mapmaking. 


of destriping mask on signal error. We plot a map of signal er¬ 
ror in 70 GHz temperature, first without a destriping mask, and 
below it with mask. The mask greatly reduces the error. 

The mask has the effect of slightly increasing the residual 
noise, as there is less data to solve the baselines from; however, 
the effect is small, well below 1 /iK, as can be seen from Fig. 
3. The effect is the opposite direction for signal error. The sig¬ 
nal error increases towards shorter baselines. The overall level 
of the signal error is still well below that of the residual noise. 
The destriping mask has a dramatic effect on the signal error, 
decreasing it typically by a factor of 2. 

We can conclude that our chosen baseline lengths are short 
enough for nearly optimal noise removal. At the same time the 
signal error is kept under control by the application of a mask. 



Fig. 6. Signal error as a function of destriping resolution for all 
LFI frequencies, for I (orange), Q (green), and U (blue). Signal 
error is computed as the rms of the difference between destriped 
and binned maps. The dashed lines in the 30 GHz case show 
results from a simulation without bandpass mismatch. 

3.2.2. Destriping resolution 

In another series of simulations we varied the destriping res¬ 
olution. The final maps were constructed at a fixed resolution 
AAg-^g=1024 in all cases. The baseline length was fixed to its nom¬ 
inal value, and a destriping mask was applied. The results are 
shown in Figs. 5 and 6. 

A higher resolution increases the residual noise slightly, 
since there are more unknowns to solve for the same amount 
of data; however, the effect is again below 1 yuK. In signal error 
the effect is in the opposite direction. In the case of temperature 
maps, increasing the resolution reduces the signal error drasti¬ 
cally, but the effect levels out near our chosen destriping resolu¬ 
tion. In polarization, the effect is less dramatic. 

The 30 GHz channel has the peculiar feature that the signal 
error is significantly higher in polarization than in temperature. 
It is also nearly completely insensitive to baseline length and de¬ 
striping resolution. The insensitivity to baseline length suggests 
that the effect is coupled to a foreground component distributed 
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30 GHz I Q U 



NW 

HUW 

NW 

HUW 

NW 

HUW 

white 

56.76 

56.78 

80.28 

80.42 

80.20 

80.41 

CRN 

15.56 

15.57 

21.75 

21.91 

21.73 

21.92 

full 

58.85 

58.87 

83.17 

83.34 

83.09 

83.32 

44 GHz 

/ 

Q 

U 


NW 

HUW 

NW 

HUW 

NW 

HUW 

white 

65.42 

65.81 

99.88 

100.44 

93.80 

94.37 

CRN 

12.52 

11.99 

19.36 

18.61 

18.28 

17.49 

full 

66.60 

66.90 

101.74 

102.15 

95.57 

95.99 

70 GHz 

/ 

Q 

U 


NW 

HUW 

NW 

HUW 

NW 

HUW 

white 

57.82 

57.95 

82.00 

82.32 

81.59 

81.94 

CRN 

3.22 

3.20 

4.50 

4.49 

4.48 

4.47 

full 

57.91 

58.04 

82.12 

82.44 

81.71 

82.06 


Table 3. Impact of detector weighting scheme on residual noise. 
We show the rms (in yuK) of binned white noise (“white”), cor¬ 
related residual noise (“CRN”), and their sum (“full”). All val¬ 
ues are averages over 10 noise realizations. NW refers to noise 
weighting, and HUW to the horn-uniform weighting scheme. 
Noise weighting minimises the total residual noise. In actual 
mapmaking we used instead horn-uniform weighting to reduce 
effects from beam shape mismatch. 


over a large fraction of the sky. The likely explanation is leak¬ 
age of synchrotron emission signal from temperature to polar¬ 
ization through bandpass mismatch. To verify this, we run yet 
another series of simulations where we artificially turned off the 
bandpass mismatch by using the same input sky model for all 
radiometers of 30 GHz. The results are shown with dashed lines 
in the uppermost panel of Fig. 6. The signal error is greatly re¬ 
duced, confirming that bandpass mismatch is the dominant fac¬ 
tor in the signal error. 

3.2.3. Horn-uniform weighting 

The purpose of the horn-uniform weighting scheme is to reduce 
leakage from temperature to polarization through systematics 
that are correlated between radiometers sharing a horn. In par¬ 
ticular this concerns leakage through beamshape mismatch. 

Changing the detector weighting changes both the binned 
and the destriped map. Therefore CRN alone is no longer a suit¬ 
able figure of merit, but we have to look at the total residual 
noise. 

Results from noise simulations are shown in Table 3. For 
completeness we show also the CRN values. We compare the 
horn-uniform weighting scheme with weighting based on the ra¬ 
diometer’s white noise level. The latter is expected to minimise 
the total residual noise. This is confirmed by our simulations. 
In all cases studied the noise weighting scheme yields a slightly 
smaller total noise residual than the horn-uniform weighting, but 
the differences are small. The white noise rms (see Table 2) typ¬ 
ically vary by 10 %, but this translates to a < 1 % effect in resid¬ 
ual noise. 

To study the effect of radiometer weighting scheme on T- 
P leakage, we performed a series of signal simulations where 
we only include the CMB temperature anisotropies in inputs. 
The resulting QIU maps consist purely of temperature leakage. 
Since foregrounds are not present, bandpass mismatch plays no 


30 GHz 

Q 

u 


NW HUW 

NW HUW 


ideal 

0.1175 

6.6e-15 

0.1385 

6.6e-15 

non-orthogonality 

0.1158 

0.0080 

0.1377 

0.0067 

realistic pointing 

0.1163 

0.0124 

0.1381 

0.0119 

realistic beam 

0.1671 

0.0890 

0.1863 

0.0893 

44 GHz 


Q 


U 


NW 

HUW 

NW 

HUW 

ideal 

0.5073 

6.5e-15 

0.5817 

6.3e-15 

non-orthogonality 

0.5239 

0.0166 

0.6016 

0.0187 

realistic pointing 

0.5313 

0.0663 

0.6081 

0.0632 

realistic beam 

0.7007 

0.2946 

0.8189 

0.3246 

70 GHz 


Q 


U 


NW 

HUW 

NW 

HUW 

ideal 

0.0731 

4.7e-15 

0.0832 

4.7e-15 

non-orthogonality 

0.0740 

0.0025 

0.0841 

0.0023 

realistic pointing 

0.0758 

0.0144 

0.0858 

0.0144 

realistic beam 

0.1900 

0.1658 

0.2177 

0.1834 


Table 4. Impact of detector weighting scheme on tempera¬ 
ture leakage. Noise weighting (NW) is compared against horn- 
uniform weighting (HUW). We show the rms (in yuK) of the Q 
and U maps from a CMB simulation with unpolarised inputs. All 
signal in 2, U maps is leakage from temperature through beam 
and pointing non-idealities. The four rows correspond to simula¬ 
tions with increasing complexity: idealised simulation (see main 
text); non-orthogonal polarization angles; realistic pointing; and 
realistic beam model. The case that most closely corresponds to 
actual mapmaking is shown in boldface. 


role, but beam shape mismatch and non-idealities in pointing are 
the only sources of leakage. We use as a simple measure of error 
in these maps the rms of the Q and U polarization maps. The 
results are shown in Table 4. 

We compare noise weighting and horn-uniform weighting in 
four simulations of increasing complexity. The first one is an ide¬ 
alized simulation, where we have eliminated any non-idealities 
that could cause leakage from temperature to polarization. We 
use identical elliptic beams and the same pointing for both ra¬ 
diometers of a horn. The polarization responses (j/^poi angle) are 
set at exactly 90° from each other, while they in reality differ 
from orthogonality by +0.5° at maximum (Planck Collaboration 
IV 2015). With these assumptions, the horn-uniform weighting 
is expected to fully cancel the temperature leakage. This is con¬ 
firmed by our simulation. 

In the second simulation we set the polarization angles to 
their real values, and in the third we add to that realistic pointing, 
which is slightly different for the M and S radiometers. Under 
these circumstances, the benefit of horn-uniform weighting as 
compared to noise weighting is still dramatic. In the last simula¬ 
tion we keep the realistic pointing and polarization angles, and 
replace the elliptic beam by the realistic main beam model. In 
this case the difference between the two weighting schemes re¬ 
duces significantly, but is still a factor of 2-3 at 30 GHz GHz and 
44 GHz. 

The level of temperature leakage reflects the beam shape 
mismatch at each frequency. Although the beams of a 44 GHz 
horn are similar in ellipticity, they have different FWHM at the 
3 % level, as contrasted to < 1 % at 30 GHz and 70 GHz (Planck 
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Collaboration IV 2015). For this reason the leakage is strongest 
at 44 GHz. 

With foregrounds present, the temperature leakage is domi¬ 
nated by bandpass mismatch. We do not attempt to correct for it 
at the mapmaking level. Bandpass mismatch is easier to model 
and account for at later stages of data processing than beam mis¬ 
match. Bandpass correction is discussed in Planck Collaboration 
II (2015). 


Our implementation of the downgrading scheme first down¬ 
grades the maps to an intermediate resolution of A^side, mid = 32. 
The Stokes I part of the map is expanded in spherical harmonics, 
treated with the smoothing beam, and the final map is then syn¬ 
thesised at the target resolution. The last downgrading step for 
Stokes Q and U maps is performed by naive averaging of higher 
resolution pixels, to minimize signal distortion. 


4. Low resolution products 

Low-resolution products are an integral part of the low-^ like¬ 
lihood. To fully exploit the information contained in the largest 
structures of the CMB sky, a full statistical description of the 
residual noise present in the maps is required. This information 
is provided in the form of pixel-pixel noise covariance matrices 
(NCVM). However, due to resource limitations they are impos¬ 
sible to employ at the native map resolution. Therefore a low- 
resolution data set is needed for the low-^ analysis. This data 
set consist of low-resolution maps and corresponding noise co- 
variance matrices. At present, the low-resolution data set can be 
efficiently used only at resolution Aside = 16 (or lower). All the 
low-resolution products are produced at this target resolution. 

We will begin by reviewing how the low-resolution maps 
were produced for the 2015 release. Then we will discuss the 
production of the 2015 LFI low-resolution noise covariance ma¬ 
trices. The noise covariance matrix section deepens the study 
performed for Planck Collaboration II (2014). 

4.1. Low-resolution maps 

Number of different approaches for producing the low-resolution 
maps exist. Previous work on production of low-resolution 
maps have been published in Keskitalo et al. (2010). An ideal 
method would simultaneously minimise the pixelization effects 
and residual noise. Realistic methods involve a trade-off between 
these two goals. 

We use the high-resolution maps described in Section 2 as 
input for the low-resolution map production. The chosen down¬ 
grading scheme is nearly identical to the one used in the previous 
2013 release (Planck Collaboration II 2014), the only difference 
being the addition of regularization noise to the final products. 

The signal in a given high-resolution pixel is weighted with 
the inverse of its noise variance when binning the signal into 
larger pixels, and hence the downgrading scheme is called the 
“noise-weighted scheme”. The weights are given by Eq. (10). 
The resulting map is equivalent to one that is directly binned 
onto the target resolution from the TOI destriped at the high res¬ 
olution. Specifically we apply to a high-resolution map the op¬ 
eration 


Pi) ' X (PhC„' Ph) /Wh = D/Wh, (20) 


where 


, _ 11 ’ if p is a subpixel of q 

(0, otherwise 


( 21 ) 


is the sum of high-resolution pixels to low-resolution pixels. 
Here subscripts “h” and “1” refer to high and low resolution ver¬ 
sion of the pointing matrix. The same matrix X also downgrades 
the pointing matrix. Pi = PhX^. After downgrading, the temper¬ 
ature component is smoothed with a Gaussian window function 
with FWHM = 440', to prevent aliasing due to high frequency 
power in the map. 


4.2. Noise covariance matrices 


The statistical description of the residual noise present in a low- 
resolution map is given in the form of a pixel-pixel noise co- 
variance matrix. For generalised destriping the formalism was 
developed in Keihanen et al. (2005) and Keskitalo et al. (2010). 
Using the formalism introduced in Sect. 2, the noise noise co- 
variance matrix N is given by 


N = 


pT(Cw + FCaF'f)'‘p 


which can be written in a dimensionally reduced form as 

N-‘ = p^c;' P + P^c;' F (F^c;* f + c-^' f^c;' p. 


( 22 ) 

(23) 


Here P, is the pointing matrix, F spreads the noise baselines to 
the TOI, Ca is the noise prior, and Cw is white noise covariance. 

For the current release, we use an implementation by 
Madam/TOAST, a Time Ordered Astrophysics Scalable Tools 
(TOAST; Kisner et al. (2015)) port of Madam. The TOAST inter¬ 
face was chosen on the basis of added flexibility. Madam/TOAST 
is capable of handling arbitrarily complex noise models. It can 
also handle bigger chunks of data, i.e., full mission versus sin¬ 
gle survey in the previously used implementation, leading to re¬ 
duced wall-clock time in computing. Some approximations and 
differences compared to high-resolution map production are in¬ 
volved, and their effects will be discussed in Sect. 4.3. These are 
the differences in baseline length and destriping resolution, as 
well as simplifying approximations. 

The above formulation in Eq. (23) describes the noise cor¬ 
relations of a map destriped and binned at the same resolution, 
according to Eqs. (6-8). For an exact description we should con¬ 
struct the matrices at resolution Aside = 1024, and subsequently 
downgrade to the target resolution. This is, however, computa¬ 
tionally unfeasible. Hence we calculate the noise covariance ma¬ 
trices at a practical intermediate resolution, which is not neces¬ 
sarily the target resolution, and subsequently downgrade to tar¬ 
get resolution. For consistency, the noise covariance matrices 
must go through the same smoothing and regularization steps 
as is applied to the low-resolution maps. 

The covariance matrix does not model the effects of using a 
destriping mask and horn-uniform weighting, which were intro¬ 
duced to reduce systematics. The impact of these approximation 
on the accuracy of the NCVM products is assessed through sim¬ 
ulations in Sect. 4.3. 

Applying equation (23) in practice will require inversion of 
a symmetric 3Apix x 3Apix matrix. The inversion is performed 
via the eigen-decomposition of the matrix. The monopole of the 
I map cannot be resolved by the map-maker, and the matrix be¬ 
comes singular. This ill-determined mode is left out of the anal¬ 
ysis. 

The chosen downgrading scheme leads to a singular covari¬ 
ance matrix, since after smoothing the number of pixels in the 
map is larger than the number of non-zero spherical harmonics, 
resulting in a significant number of zero eigenvalues. We regu¬ 
larise the products by adding a small amount of white noise to 
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Fig. 7. Reduced statistics for simulated full mission 30 GHz 
data. Initial resolution of the noise covariance matrix is A^side = 
32. The level of idealization decreases from top to bottom. 


the maps, and a corresponding diagonal covariance to the matri¬ 
ces. The level of regularization noise is chosen to be 2 juK RMS 
for /, and 0.02 juK RMS for Q and U at A^side = 16 resolution. 

As seen from equation (23), the noise covariance compu¬ 
tation requires two inputs: the detector pointings and the noise 
model derived from the flight data. 

For the noise covariance matrix production we use the most 
representative noise models available, namely the FFP8 noise es¬ 
timates (Planck Collaboration XII2015). The noise model com¬ 
prises of daily 1// model parameters. 

The NCVM production uses the FFP8 pointing solution of 
TOAST as an input. It reproduces the data processing centre 
(DPC) pointing to sub arc second accuracy except for a few iso¬ 
lated time periods where the flagging in TOAST differs from the 
DPC one. The difference between DPC and FFP8 bitmaps for all 
LFI frequency channels are illustrated in Fig. 9, and the statistics 
of hits map differences are summarised in Table 5. For survey- 



0.90 0.95 1.00 1.05 1.10 1.15 1.20 


Reduced 

Fig. 8. Reduced statistics for simulated full mission 30 GHz 
data. Initial resolution of the noise covariance matrix is varied 
from Aside = 16 to 64. 


by-survey comparisons see Planck Collaboration XII (2015); 
this level of discrepancy is negligible for the results. The dif¬ 
ference in hit counts is greatest around the ecliptic poles, where 
we have the most hits per pixel. 

For noise covariance computations the output of 
Madam/TOAST consists of inverse noise covariance matri¬ 
ces defined by Eq. (23), specifically one inverse matrix per 
radiometer for a given time period. For all LFI frequency 
channels we use 0.25 s baselines in the matrix production; the 
exact numbers are given in Table 5. The matrices are computed 
at the highest feasible resolution of Aside = 64. The individual 
inverse matrices are merged together to form the actual inverse 
NCVMs. These inverse matrices are then inverted using the 
eigen-decomposition of a matrix, and subsequently downgraded 
to the target-resolution using the same downgrading scheme as 
applied for the maps. The 70 GHz full mission noise covariance 
matrix computation took approximately 6 hours of wall-clock 
time and 23 000 CPU hours. 

4.3. Validation 

We have performed a series of simulations to assess the effect of 
different parameters and approximations to the quality of noise 
covariance matrices. 

We performed noise MC simulations with the FFP8 pipeline 
(Planck Collaboration XII 2015), but by varying the mapmak¬ 
ing parameters. For all combinations of simulation parameters 
we have generated 100 realizations of noise-only maps. We al¬ 
tered the baseline length, the destriping resolution, the weight¬ 
ing scheme, or the inclusion of destriping mask. In the noise 
covariance matrix production we altered the baseline length and 
destriping resolution. We chose to perform the simulations for 
the 30 GHz channel, which requires the least computational re¬ 
sources. Figures 7 and 8 summarise the results from our simula¬ 
tions. 

We calculate for each simulated noise map m the reduced 
statistics, defined as 




2 


^dof 


(24) 
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Table 5. Frequency-specific information related to the noise covariance matrix production. From left to right: the statistics of 
difference between DPC and FFP8 full mission bitmaps for all LFI frequency channels; the statistics of DPC full mission hits maps; 
and the baseline lengths used in the noise covariance matrix production, both in seconds and as a number of samples. 


Channel 

Difference hitmap statistics 

Hitmap statistics 

Baseline length 

Min 

Max 

St. Dev. 

Mean 

Min 

Max 

[S] 

Samples 

30 GHz 

-18 

18 

0.71 

1180 

172 

92100 

0.2461 

8 

44 GHz 

-22 

28 

1.02 

2530 

574 

134000 

0.2578 

12 

70 GHz 

-34 

44 

1.88 

8560 

2130 

266000 

0.2539 

20 




-5 Hits 5 

Fig. 9. The difference of DPC and FFP8 full mission bitmaps 
for all LFI frequency channels. The maps are shown in galactic 
coordinates. Top: 30 GHz. Middle: 44 GHz. Bottom: 70 GHz. 


where N is the noise covariance matrix, and the number of de¬ 
grees of freedom A^^of = 3A/pix. In the presence of ill-conditioned 
eigenmodes, the degrees of freedom are reduced accordingly. 
For unregularized smoothed noise covariance matrices the ef¬ 
fective number of degrees of freedom equals the number of non¬ 
zero eigenvalues. The distribution of reduced will peak at 1 
for a noise covariance matrix that models perfectly the properties 
of the simulated noise maps. 

In Fig. 7 the degree of idealization reduces from top to bot¬ 
tom. We add one effect at a time, so that the last panel repre¬ 


sents a fully realistic simulation. We start from an idealized case 
(panel A), where both the noise maps and noise covariance ma¬ 
trices have been calculated at resolution Aside = 32, with one 
sample baselines (0.03125 s). In panel B we downgrade both 
maps and matrix to the resolution of Aside = 16, and in panel C 
we regularize the maps and matrices. We observe that the effects 
of downgrading and regularization on the;^^ result are marginal. 

The key parameter in noise covariance matrix computation is 
the baseline length. In earlier studies we have demonstrated that 
shorter baseline length in the noise covariance matrix production 
models the residual noise better (Planck Collaboration II 2014). 
Given the enhanced computing resources it was possible to bring 
the baseline length to one sample. Panels D - F summarize the 
results. In panel D we increase the baseline length in MC sim¬ 
ulations from one sample to 8 samples (0.25 s), which was the 
value used in actual mapmaking. Panels E-F show the effect of 
increasing the baseline length in matrix computation from one 
sample to first 2 samples, and then to 8 samples (0.25 s). 

We find that reducing the baseline length in NCVM compu¬ 
tation below 0.25 sec does indeed improve the results, but at the 
same time the computer memory requirements increase rapidly. 
For the 2015 release we have used 0.25 sec baselines for all LFI 
frequency channels, but we propose the adoption of even shorter 
baselines in future releases, resources permitting. 

Panel G shows the effect of changing the destriping resolu¬ 
tion in MC simulations from Aside = 32 to the realistic value 
Aside = 1024. 

The noise covariance computation makes two further devi¬ 
ations from the high-resolution mapmaking: the horns are not 
uniformly weighted (panel H) and it does not take into account 
the destriping mask (panel I). Panel I finally represents a fully 
realistic simulation. The effects in panels G-I are shown to be 
much smaller than the effects of the baseline length and the de¬ 
striping resolution. 

We examine the impact of resolution further in Fig. 8. As 
motivated in previous section the matrices should be calculated 
at the highest computationally feasible resolution. We calculated 
a few noise covariance matrices with varying initial resolution. 
Each matrix was downgraded to resolution Aside = 16 and com¬ 
pared against a downgraded set of FFP8 noise maps (100 re¬ 
alisations). As expected. Aside = 64 gives the best agreement. 
Increasing the initial resolution beyond Aside = 64 is likely to 
improve results further, but the expected improvement would be 
a very small for a huge computational effort. Increasing the ini¬ 
tial resolution beyond Aside = 64 to Aside = 128 will increase the 
matrix size by a factor of 16 to 2.5 TB. 
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Fig. 10. LFI 30 GHz channel maps. From top to bottom: tem¬ 
perature; 2, and U polarization. The temperature map has res¬ 
olution A^nside=1024. The polarization maps have been down¬ 
graded to resolution A^side=256 and smoothed with a FWHM=1° 
Gaussian beam. The polarization maps are not corrected for 
bandpass mismatch leakage. 

5. Results 

5.1. High-resolution maps 

The 2015 release includes a large number of map products. A 
complete list is given Planck Collaboration II (2015). Figures 
10-12 show the main products, which are the frequency maps 
for the complete four-year data set. Shown in these figures are 
the temperature map, and Q and U polarization. For illustration 
purposes the polarization maps have been smoothed with a sym¬ 
metric Gaussian beam with FWHM=1°. 

The released maps are in Galactic coordinates. The division 
of the polarization signal into Q and U components is depen- 


-20 f^K 20 

Fig. 11. LFI 44 GHz channel maps. Fom top to bottom: temper¬ 
ature, Q, and U polarization. 


dent on the chosen coordinate frame. The Q map represents a 
component where the polarization direction is aligned with the 
local meridian (positive Q) or with the local latitude (negative 
Q). The U map represents a component where the polarization 
direction is at a 45° angle with respect to the meridian. The polar¬ 
ization maps shown here are subject to bandpass leakage (Planck 
Collaboration II 2015), which must be taken into account in the 
subsequent analysis. 

We do not correct for beam shape in the mapmaking process. 
The shape of a point source in a map directly refiects the beam 
shape. In particular, the power that is lost into the far sidelobes 
(Planck Collaboration IV 2015), is missing from the point source 
as well, and must be taken into account when determining source 
fiuxes. 

As auxiliary information we provide for each high-frequency 
map an estimate of the white noise covariance. The matrix gives 
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Fig. 12. LFI 70 GHz channel maps. From top to bottom: temper¬ 
ature, 2, and 1] polarization. 


the noise correlation between Stokes components within each 
pixel. We provide also a corresponding hit count map for each 
map. 

5.2. Sky coverage 

Among other map products are single-survey maps that include 
data from one 6 month survey, and single-year maps than in¬ 
clude two surveys each. These do not give full sky coverage. As 
discussed in Sect. 2.1.3, construction of (/, Q, U) Stokes com¬ 
ponents requires that a pixel is scanned by a pair of horns, with 
complementary polarization sensitivity. Pixels that are scanned 
by a single horn or radiometer, can be used for destriping, but 
not for the final map. 

The sky coverage fraction of the frequency maps are given 
in Table 6. We present three numbers for each map. In the first 


column (“map”) we give the sky coverage of the final map. This 
is the fraction of sky that has been scanned by two radiometer 
pairs so that we have been able to solve the full /, Q, U triplet. 

The second column (“scan”) gives the fraction of sky that 
is scanned by at least one radiometer, and is thus available for 
the solution of the destriping equation. The focal plane is de¬ 
signed in such a way that in most cases a pair of horns follows 
the same scanning path. Each pixel is thus covered by four ra¬ 
diometers, and the difference between the two columns is small. 
An exception to this is the 44 GHz frequency channel, where 
horn LFI-24 does not have a counterpart. A significant fraction 
of the sky in single surveys is covered by horn LFI-24 alone. As 
a consequence, the sky fraction available for destriping is signif¬ 
icantly higher than that available for the final map. In the third 
column (“destriping”) we give the sky fraction that remains after 
the galactic mask is applied. This is the fraction eventually used 
for destriping. 

Single-survey temperature maps for each full frequency are 
shown in Fig. 13. We indicate by different colour scales the re¬ 
gions used for the final map or for destriping. The combination 
of blue and purple regions is included in the final map. The size 
of this region is given in the “map” column of Table 6. The 
Galactic region, shown in blue, is masked out in the destriping 
phase, and not used for solving the baselines. The red colour, 
mostly visible at 44 GHz, indicates the region that is scanned by 
a single horn, and is used for destriping, but is not included in 
the final map. The purple and red regions together correspond to 
the column headed “destriping” in Table 6. The region shown in 
white does not contain useful data. 


5.3. Low-resolution products 

The 2015 release includes low-resolution maps and matrices 
both for the full mission and for a data selection covering sur¬ 
veys 1,3, and 5-8 (Planck Collaboration II 2015) at resolution 
of ^side = 16. They have been processed with the noise-weighted 
downgrading scheme, and subsequently the Stokes I component 
has been smoothed with a Gaussian beam of FWHM = 440'. 
Finally the products have been regularised by adding 2 juK RMS 
for /, and 0.02 pK RMS for Q and U of white noise at Aside = 16 
resolution. 

Figure 14 illustrates the full mission low-resolution Stokes 
/, Q and U maps for each Planck LFI frequency channel. Figure 
15 illustrates a single column of the 70 GHz full mission noise 
covariance matrix. The column corresponds to a pixel number 
1248 in the HEALPix nested pixelization scheme for Aside = 16 
resolution. The pixel values in the plot represent correlation co¬ 
efficients, since each pixel value has been normalised by 

^{m^){mp. In this normalization the reference pixel automati¬ 
cally gets unit value, and is later set to zero to bring out finer 
details of the noise covariance matrix. 


5.4. Characterization of noise 

5.4.1. Residual noise at high multipoles 

Residual noise at high multipoles is dominated by non-uniform 
white noise, which is uncorrelated between pixels. We provide 
a description of the white noise component in the form of the 
white noise covariance, computed according to Eq. (14). The co- 
variance matrix consists of a symmetric 3x3 matrix for every 
map pixel. The matrix gives the noise correlation between Stokes 
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Fig. 13. Sky-coverage of single-survey temperature maps for 30, 44, and 70 GHz. Left: 30 GHz. Middle: 44 GHz. Right: 70 GHz. 
From top to bottom. Surveys 1 to 8. The combination of red and purple regions is used for destriping. The combination on purple 
and blue regions is included in the final map. Unobserved regions are shown in white. 
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Table 6. Sky coverage (%) of full, single-year, and single-survey frequency maps. We give the fraction of the sky covered by the 
final map (“map”), the fraction that was scanned by at least one horn (“scan”) and is available for destriping, and the fraction that 
remains after applying a galactic mask and that was used for destriping (“destriping”). The coverage of the destriping mask can be 
read from the first row of the “destriping” column.) 


Survey 

PID range 


30 GHz 



44 GHz 


70 GHz 

map 

scan 

destriping 

map 

scan 

destriping 

map 

scan 

destriping 

full 

3-44072 

100.0 

100.0 

78.74 

100.0 

100.0 

89.37 

100.0 

100.0 

89.67 

year 1 

3-10957 

99.999 

100.0 

78.74 

99.998 

100.0 

89.37 

100.0 

100.0 

89.67 

year 2 

10958-21482 

99.98 

100.0 

78.74 

99.998 

100.0 

89.37 

100.0 

100.0 

89.67 

year 3 

21483-32761 

99.72 

99.74 

78.66 

99.91 

100.0 

89.37 

99.65 

99.65 

89.38 

year 4 

32762-44072 

95.68 

95.76 

76.61 

96.87 

98.68 

88.35 

97.38 

97.43 

87.63 

survey 1 

3-5483 

97.20 

98.11 

76.88 

93.93 

98.03 

87.41 

97.94 

97.95 

87.62 

survey 2 

5484-10957 

97.48 

97.58 

77.40 

93.31 

98.13 

87.89 

97.47 

97.49 

87.62 

survey 3 

10958-16454 

97.62 

97.71 

76.49 

93.65 

98.23 

87.61 

97.61 

97.63 

87.30 

survey 4 

16455-21482 

91.88 

91.98 

72.37 

89.53 

94.92 

84.85 

92.40 

92.43 

82.89 

survey 5 

21483-27404 

90.89 

91.02 

70.27 

88.43 

96.70 

86.08 

92.44 

92.50 

82.33 

survey 6 

27405-32761 

87.79 

88.06 

72.35 

86.10 

94.56 

85.29 

89.95 

90.18 

82.04 

survey 7 

32762-38574 

85.40 

85.60 

68.29 

83.70 

92.49 

82.80 

88.43 

88.52 

79.65 

survey 8 

38575-44072 

80.01 

80.26 

67.54 

78.92 

89.61 

81.15 

83.83 

84.09 

77.05 



Fig. 14. LFI full mission low resolution maps, at A^side = 16. From left to right: intensity /; polarization Q component; and polariza¬ 
tion 1] component. Top\ 30 GHz. Middle: 44 GHz Bottom: 70 GHz. Units are in //Kcmb- 


components, /, 2, and U, within each pixel. White noise gives a 
lower limit for the full residual noise level. 

Table 7 shows the average white noise covariance per 
pixel in the frequency maps. We show the average value of 
II, IQ, lU, QQ, QU,UU noise covariance, as average over all 
pixels. We also give the expected white noise level in the angular 
power spectrum, calculated as 

Cn = ( 25 ) 

2Vpix 


and 

CTe = C7t = l((QQ) + (26) 

where A^pix is the number of pixels in the map (12 582 912) and 
brackets indicate an average over pixels. 

The 44 GHz channel shows an imbalance between Q and U 
components. This is due to horn LFI-24, which does not have 
a counterpart with complementary polarization sensitivity. The 
same effect can be seen in the CRN plots in Sect. 3. 
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Fig. 15. A single column of the 70 GHz full mission noise covariance matrix. The reference pixel is pixel 1248 in the HEALPix 
nested pixelization scheme. The noise covariance matrix has been normalised to show correlation coefficients. The reference pixel 
is set to zero to highlight more of the finer details of the correlation structure. 


Table 7. White noise level for frequency maps, computed from the white noise covariance matrix. We show the average white noise 
covariance per pixel for Stokes components, and white noise spectrum level derived from the same numbers. All numbers are in 
units. 


Channel 

II 

QQ 

UU 

IQ 

lU 

QU 

Ctt 

Cee^ Cbb 

30 GHz 

3224.4 

6467.8 

6465.1 

-19.2 

-71.0 

-41.2 

0.003220 

0.006458 

44 GHz 

4331.0 

10088.0 

8906.9 

-2.4 

24.3 

-136.9 

0.004325 

0.009485 

70 GHz 

3358.5 

6775.7 

6713.7 

28.1 

-69.0 

-83.3 

0.003354 

0.006736 


We use half-ring maps to derive an estimate for the level of 
correlated residual noise at high multipoles. We find that even af¬ 
ter destriping there remains correlated residual noise at the high¬ 
est multipoles, and it affects the auto-spectra at the 1% level. 
In Fig. 16 we plot the estimated noise spectrum in the high- 
multipole regime, as a fraction of the white noise level reported 
in Table 7. We show the mean of the spectrum over multipoles 
^ = 1150 - 1800. We show in the same figure also the mean 
of 10 000 FFP8 noise Monte Carlo realizations, all derived from 
same noise spectrum. The error bars on the MC estimated indi¬ 
cate the statistical Icr variation of the realizations, not the uncer¬ 
tainty of the noise model. 


5.4.2. Residual noise at low multipoles 

We have three complementary methods for chracterization of 
correlated residual noise at low multipoles: the noise covariance 
matrix; MC simulations; and the half-ring noise map. We now 
compare the noise estimates from the three methods against each 
other. 

Figure 17 shows the reduced statistics for full mission 
noise covariance matrices. The x^ tests have been carried out 
with 10 000 FFP8 noise realizations with realistic mapmaking 
parameters. The vertical lines indicate the reduced value for 
the half-ring noise maps. We observe that at the 70 GHz fre¬ 
quency channel, which is the most important for the cosmologi¬ 
cal parameter estimation, the noise covariance matrix models the 
residual noise the best. 
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Frequency channel 


Fig. 16. Comparison of noise estimates at high multipoles, as a 
fraction white noise level (blue). The green dots show the esti¬ 
mate from half-ring maps. The orange dots show the mean of 
10000 Monte Carlo simulations. The error bars include the sta¬ 
tistical variation only. The values shown are mean values of the 
noise spectrum in the range i = 1150-1800. The absolute white 
noise levels in are reported in the last two columns of Table 
7. 
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Fig. 17. Reducedstatistics for 2015 full mission noise covari¬ 
ance matrices, t^sts have been carried out with 10000 FFP8 
noise realizations. Vertical lines mark the reduced value for 
half-ring noise maps. 


Figures 18-20 show noise bias estimates for full frequency 
maps at low multipoles. We show the six noise spectra (TT, 
EE, BB, TE,TB, EB) constructed from the NCVMs, and corre¬ 
sponding estimates from FFP8 noise Monte Carlo simulations. 
We show the average and median over 10 000 noise realizations. 
We plot in the same figure also the spectrum of the half-ring 
noise map. As this represents one realization only, it has much 
stronger variations than the other estimates. 


6. Conclusions 

In this paper we have described the mapmaking procedure used 
for the 2015 data release from Planck LFI, and characterised the 
noise in the map products. For the first time, we have shown 
polarization Q and U maps at LFI frequencies. 

High-resolution maps were produced with the Madam map¬ 
making code. The algorithm is based on a destriping technique, 
which is extended to short baseline lengths with the help of a 
noise prior. We validated the procedure through noise and signal 
simulations. We studied the effect of two main parameters that 
control the mapmaking procedure, namely the baseline length, 
and the destriping resolution. We showed that the values chosen 
for actual pipeline are appropriate for the data set and allow for 
good noise removal, while keeping the signal error at low level. 

We have paid special attention to the reduction of unwanted 
systematic effects. We applied a Galactic mask to reduce signal 
error arising from bandpass mismatch and strong signal gradi¬ 
ents. We assessed the impact through simulations, and showed 
that the signal error is indeed reduced significantly, while the 
cost in residual noise is low. 

The horn-uniform weighting scheme was introduced to re¬ 
duce leakage of temperature signal to polarization through 
beamshape mismatch. The effect at 30 GHz and 44 GHz was 
significant, less so at 70 GHz. In all cases, however, the cost in 
residual noise level was small. 

The polarization map products presented here are subject to 
bandpass leakage. The leakage correction applied is discussed in 
Planck Collaboration II (2015). 

The maps are not corrected for beam shape, or for the power 
to the far sidelobe part of the beam. In particular, images of point 
sources reflect the radiometer beam shapes. The exact beam 
shapes must be accounted in the estimation of source fluxes from 
LFI maps. 

We have also described the production of low-resolution 
maps and noise covariance matrices. The low-resolution maps 
are produced by downgrading the high-resolution maps to the 
target resolution of Aside = 16 with the noise-weighted scheme, 
and the I component is subsequently smoothed with a Gaussian 
beam. The noise covariance matrices are initially calculated at 
the resolution of Aside = 64 with Madam/TOAST software, and 
subsequently the same downgrading scheme is applied to them 
as for the maps. Finally both the low-resolution maps and matri¬ 
ces are regularised by adding a small amount of white noise to 
all products. 

The matrices were validated against noise Monte Carlo sim¬ 
ulations. The match between the matrices and the simulations 
was best at 70 GHz, which is the most important channel for cos¬ 
mology analysis. We performed a series of simulations to quan¬ 
tify the effect of each approximation involved in the noise co- 
variance computation. The most important effects turned out to 
be the high destriping resolution, and the limited baseline length. 

In addition to the noise covariance matrices, we charac¬ 
terised the residual noise in maps through noise Monte Carlo 
simulations and the noise maps constructed from half-ring maps. 
The noise estimates from the complementary methods show 
good agreement. 

Acknowledgements. The Planck Collaboration acknowledges the support of: 
ESA; CNES, and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF 
(Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MINECO, 
JA and RES (Spain); Tekes, AoF, and CSC (Finland); DLR and MPG 
(Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); 
RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); ERC and PRACE 
(EU). A description of the Planck Collaboration and a list of its mem¬ 
bers, indicating which technical or scientific activities they have been in- 


19 


























































Planck Collaboration: LFI mapmaking 



Multipole ^ 

Fig. 18. Noise bias levels at 30 GHz. Noise bias from 2015 full mission noise covariance matrix is plotted in orange, the noise MC 
statistics from 10000 FFP8 simulations are plotted in green (for the mean), blue (median) and grey (+1 sigma region), as well as 
the noise estimate of the half ring noise map in purple. 
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Fig. 19. Noise bias levels at 44 GHz. 
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Fig. 20. Noise bias levels at 70 GHz. 


volved in, can be found at http://www.cosmos.esa.int/web/planck/ 
planck-collaboration. Some of the results of this research have been 
achieved using the PRACE-3IP project (FP7 RI-312763) resource Sisu based in 
Finland at CSC. We thank CSC - IT Center for Science Ltd (Finland) for compu¬ 
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DE-AC02-05CHI 1231. This work has made use of the Planck satellite simu¬ 
lation package (Level-S), which is assembled by the Max Planck Institute for 
Astrophysics Planck Analysis Centre (MPAC) Reinecke et al. (2006). Some of 
the results in this paper have been derived using the HEALPix package Gorski 
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